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ABSTRACT 

The dependence of microlensing time scale frequency distributions and optical depth 
toward the galactic center on galactic model parameters is explored in detail for a distri- 
bution of stars consisting of the Zhao (1996) bar and nucleus and the Bahcall and Soneira 
(1980) double exponential disk. The high sensitivity of these two microlensing measures 
to the circular velocity model, velocity dispersions, bulge mass, direction of the line of 
sight, bar axis orientation, star spatial distribution and the stellar mass function means 
no single galaxy property can be constrained very well without constraining most of the 
others. However, this same sensitivity will make microlensing a powerful member of the 
suite of observational techniques that will eventually define the galaxy properties. The 
model time scale frequency distributions are compared with that determined empirically 
by the MACHO group throughout. Although the MACHO empirical data are matched 
quite well with a nominal velocity model and with a mass function only of hydrogen 
burning stars that varies as m~^'^ to m~^'^ in the M star region, uncertainties in galactic 
structure, kinematics and content together with the paucity of published microlensing 
data preclude any claim of the model representing the real world. A variation of the 
mass function ~ m^^ in the M star region obtained from recent star counts, both local 
and in the galactic bulge, fails to yield a sufficient number of short time scale events 
compared to the MACHO data. The high sensitivity of the microlensing measures to 
the direction of the line of sight may mean that sufficient microlensing data is already 
in hand to constrain the bar distribution of stars. The procedure developed here for 
determining the time scale frequency distribution is particularly convenient for rapidly 
incorporating model changes as data from all sources continues to accumulate. 



1. Introduction 

The successful development of techniques to mon- 
itor microlcnsing events toward the center of the 
galaxy has provided an additional observational probe 
of galactic properties. Already an optical depth that 
exceeds predictions (Alcock et al 1996; Udalski et al 

1994) has confirmed other indications (Dwek et al 

1995) that the galactic bulge must be bar shaped 
with the long axis pointing toward us (Paczynski et al 
1994). More massive lenses tend to have longer time 
scales, so microlensing places statistical constraints on 
the mass function of the stellar lenses, although other 
variables entering into the time scale weaken these 
constraints. Attempts to define the mass function of 
the lenses from statistical arguments or from the event 
time scale frequency distribution have led to models 
with a large population of brown dwarfs in the bulge 
(Han and Gould 1996; Han 1998) but other models 
exist with few or no brown dwarfs (Zhao, Rich & 
Spergel 1996; Zhao and deZeeuw 1997; Mera, Cabrier 
& Schaeffer 1998). By measuring optical depth and 
the distribution of event time scales as a function of 
the line of sight direction, microlensing can also con- 
strain the distribution of stars in the bar. So far there 
is insufScient data published [40 secure events from 
the MACHO (MAssive Compact Halo Objects) group 
(Alcock et al 1996) and 12 from the OGLE (Optical 
Gravitational Lensing Experiment)] group (Udalski et 
al 1994)] for such constraints to be very definitive. In 
fact, the lack of sufficient observations of all types 
that can constrain the properties of the galaxy means 
that few of the potential constraints on galactic struc- 
ture, kinematics or content that microlensing can con- 
tribute to can be yet realized — the major exception 
perhaps being the confirmation of the bar. 

The microlensing optical depth is the probability 
that a microlensing event will be occurring at any in- 
stant in a particular direction for a single source. The 
time scale frequency is the number of events of a par- 
ticular duration that occur per unit time, where the 
unit of time can be the entire observational period 
and the number of sources can be all of those mon- 
itored. These two measures are routinely obtained 
from a microlensing observational program aimed at 
monitoring stars toward the galactic center, and form 
the basis for constraining galactic properties. It is 
important to appreciate the sensitivity of such data 
to the various galactic parameters on which they de- 
pend, since they are related to the galaxy through 



models with assumed vahies of the parameters. Most 
of the galactic parameters such as stellar velocities 
as a function of position are poorly constrained ob- 
servationally, which may allow several ways to match 
microlensing data with very different galactic models. 
Given this uncertainty, we wish to investigate the sen- 
sitivity of time scale frequency distributions and op- 
tical depths to changes in the values of several galac- 
tic parameters, and compare the model microlensing 
measures at every step with the empirical measures 
deduced from the microlensing events followed by the 
MACHO group during the summer of 1993 (Alcock 
et al 1996). We eliminate the two stars that are likely 
to be variable stars, two stars that might fail to pass 
the MACHO criteria for inclusion and one binary star 
from the data set. This leaves a rather meager data 
set of only 40 events, but we chose not to add the 12 
OGLE events in the interest of having a data set all of 
whose entries are analyzed the same way. The small 
gain in the statistics from adding the OGLE events 
does not warrant the added complication. 

The best model of the distribution of stars in the 
galactic bulge is that due to Zhao (1996). This model 
of the bar-shaped galactic bulge is based on a model 
of Dwek et al (1994) that reproduces the CORE IR 
surface brightness distributions, but which is also con- 
sistent with all other available observations of bulge 
kinematics. A highly peaked galactic nucleus is added 
to the triaxial "boxy" Gaussian distribution of the bar 
to produce the observed central kinematics. We use 
the Zhao distribution of stellar mass (bar plus nu- 
cleus) and the double exponential disk distribution 
of stars of Bahcall and Soneira (1980) as a basis for 
determining microlensing event time scale frequency 
distributions and microlensing optical depths as func- 
tions of galactic parameters and direction of the line 
of sight. 

In Section 2 we develop a general scheme for de- 
termining time scale frequency distributions, in which 
it is relatively easy to substitute a variety of stellar 
spatial distributions (although limit thc^se partic- 
ular changes to different bar orientations), stellar ve- 
locities that depend on position in the galaxy, mass 
functions, luminosity functions and line-of-sight di- 
rections. In Section 3 we discuss galactic models in 
general and give details of the particular model that 
will be the basis for the calculations. A short Sec- 
tion 4 describes the Monte Carlo integration scheme 
for two choices of independent variables in the six di- 
mensional integral defining the time scale frequency 



2 



distribution. Section 5 is a scries of example calcula- 
tions, where the first example is a demonstration of 
the robustness of the time scale frequency distribution 
obtained with either choice of independent variables. 
The remaining examples show the effects of changing 
the visibility of the sources, the orientation of the long 
axis of the galactic bar relative to the line from the 
observer to the galactic center, the velocity model, the 
velocity dispersions within one model, the mass func- 
tion and the direction of the line of sight. The empiri- 
cal time scale frequency distribution is superposed on 
most of the graphs for comparison. The last example 
in Section 5 demonstrates an almost perfect match of 
the average of model time scale frequency distribu- 
tions over the 24 MACHO field directions with the 
empirical distribution. 

A brief derivation of the optical depth averaged 
over the source distribution along the line of sight be- 
gins Section 6. This expression is used to demonstrate 
a relatively large variation in optical depth over a re- 
gion that includes Baade's window. The spatial dis- 
tribution of optical depths along each MACHO line of 
sight is compared with the average optical depth de- 
termined empirically from the MACHO observations, 
where the same model giving such a good match to 
the empirical time scale frequency distribution gives 
an average optical depth that is about 2cr below the 
empirical mean. The average of the model optical 
depths over the MACHO fields is raised almost to the 
empirical mean if the bulge mass is multiplied by 1.5. 
The variation of optical depth for a particular direc- 
tion on source visibility and orientation of the bulge 
axis is also determined. In both Sections 5 and 6, 
we show the contributions of each separate part of 
the distribution of stars, which demonstrates how lit- 
tle the disk contributes to both the event rate and 
optical depth relative to the bulge. The discussion 
in Section 7 points out virtues of and weaknesses in 
our particular formulation, emphasizes the sensitivity 
of results to assumptions within the models, cautions 
the reader about taking our model match to the data 
too seriously given the large uncertainties in the real 
galactic model, but points out the eventual power of 
microlensing to play a major role in the entire suite 
of observations that will define the properties of the 
galaxy in the future. Our conclusions are enumerated 
in a brief Section 8. 



2. Event time scale frequency distribution 

A microlensing event is said to occur whenever a 

source star and lens star pass each other at an angular 
separation within that of the Einstein ring radius Re 
of the lens. The time scale for such an event is defined 
as Ie = Re/v, where v is the magnitude of the relative 
transverse velocity between source and lens projected 
onto the lens plane. This projected transverse relative 
velocity is given by 



V = - Vo - (vs - Vo 



Dr 



Dr 



(1) 



where Vl, Vg and Vq are the transverse velocities of 
the lens, source and observer respectively relative to 

inertial space, and Dol and Dos arc the distances to 
the lens and source. The Einstein ring radius is 
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where G is the gravitational constant, m is the lens 
mass and c is the velocity of light. 

A single source at distance Dqs will lead to mi- 
crolensing events with lenses within mass range dm 
of mass m, within dDoL of Dql, and with Einstein 
ring radii Re at the rate of 
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where {dni,/dm)dm is the number density of lenses 
with masses within dm of m so that {driL /dm,)dDoLdm 
is the areal density of lenses in the mass range dm in 
the slab of thickness dDoL at Dql- The rate is ex- 
pressed in terms of mass density of lenses through 
dpz^{m,DoL)/dm, = mdn,^(m, DoL)/dm, where the 
m in the denominator cancels the m in i?^ in the 
second form. The transverse velocity v has compo- 
nents and Vf in the directions of increasing galactic 
latitude and longitude respectively, such that 



(vl + vl)'/' 



(4) 



with 
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'i^LSR + 
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In Eqs. (||), Vlsr is the circular velocity of the local 
standard of rest, Vob,e are the peculiar velocity com- 
ponents of the observer relative to the LSR, Vlo{Dol) 
is the circular velocity of the galaxy at Dql, v'^^ ^ are 
the random lens velocities relative to the mean galac- 
tic rotation, f so (-Dos) is the circular velocity of the 
galaxy at fosi and -Wg^^ are the random source veloc- 
ities. We shall separate the distribution of stars in the 
galaxy into disk and triaxial bulge distributions with 
their own circular velocity distributions and velocity 
dispersions, where the latter are Gaussian distributed 
about the mean circular velocities. 

The number of source stars in a given field that 
are visible in range dDos at Dqs is proportional to 
ns{Dos)D'^^^ , where ns{Dos) is the number density 
of sources at Dqs, and where the exponent follows 
from the increase in area of the cross section of the 
angular field with distance and from the assumption 
that the fraction of stars with luminosities greater 
than some varies as (Kiraga and Paczyhski 
1994). Hence, the fraction of the visible sources in 
the field that are within the slab of thickness dDos is 
given by 



n{Dos)Dl+''^dDos/ n{Dc 



)D 
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dDr 



The expression is constrained to lenses within 
dDoL at Dql and to lens masses within dm, of m. To 
further constrain the event rate to particular ranges 
of lens and source random velocities and to a range 
of source distances dDos, we must multiply by the 
respective fractions of lenses and sources contained 
within the small range of each parameter. Hence 



dF 
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mass within dm of m. The Gaussian velocity distri- 
butions are represented by 



1 



■ exp 



where Oxy is the rms value of v'^^ The circular ve- 
locities of the galaxy are assumed known from the 
galactic model. The particular values of all seven pa- 
rameters in Eq. determine a particular value of 
the time scale Ie = Re/v which we have substituted. 

To determine the total event rate per unit time 
scale interval for a particular time scale Ie = t^, we 
multiply Eq. by Sits — i'e) and integrate over 
the seven dimensional volume. Before doing this, it 
is expedient to normalize all distances by I?8 = 8 
kpc, the distance to the galactic center, all velocity 
components by v^sr, and the masses by the solar mass 
Mq. The time scale Ie is normalized by 



to 



IAGAU.,D^ 



66.72 days 



(7) 



which is the time scale for an event with a solar mass 
lens located at a distance of 4 kpc with the source at 
the galactic center (8 kpc). The source is assumed 
at rest and the lens velocity and observer velocity are 
both assumed to be the circular velocity v^sr consis- 
tent with a flat galaxy rotation curve. The numerical 



value for to is obtained with v. 



210 km/sec (Al- 



varez, May & Bronfmann 1990; Clemens 1985). These 
normalizations yield 
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as the number of events/sec/source/ (unit Ie) at time 
scale Ie = t'^. The numerical coefficient means is 



j n{Dos)Dl+'^'^dDos (6) expressed in M^/pc^, and C = Dos/D„ z = Doi^jD^, 
J and it is understood that all velocities are normalized 



is the rate of events for a single source within dDos of 
Dqs that also has random velocities within dug^dwgj 
of (v'gf^iV'sf) for lenses within dDoi, of I^ol with ran- 
dom velocities within dv'j^^dv'j^g of (w^^,w^f) and with 



and it is understood that all velocities are normalized 
by Vlsr and Ie by tg. The integral limits correspond 
to the variables in the same order as the differen- 
tials in the integrand. The integral over starts at 
0.1 instead of to avoid the singularity in the inte- 
grand. There are likely to be no sources involved in 
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microlensing events closer to us than 0.8 kpc. The up- 
per limit on ^ extends the range beyond the galactic 
center were source visibility begins to fade from ei- 
ther distance or blending. One can experiment with 
the effect of changing these limits, since the extent of 
blending of stellar images depends on observing site 
seeing. 

The integral is to be evaluated approximately with 
a Monte Carlo technique, but first the dimension is re- 
duced by one by changing variables from m, z, ^, w^^, 

'"'si,' '"'si to^B' 2:, C, '"'li' '"'sb' ^'si and integrat- 
ing over the 5 function. The Jacobian determinant of 
the transformation is simply dmldt^- In dimension- 
less units, 



such that 



dm 

dt: " 



/ mz{(^ — z) 



(9) 



(10) 



z{C-z) 

The change of variables in the integrand of Eq. (^) 
is effected by replacing dm with {dm/dtE)dtE. Inte- 
gration over fixes at the value selected by the 5 
function, and the remaining six dimensional integral 
is 



16 X 1033 



I GDI 



S^LSR 



(3.0856 X 1018)3 

00 poo poo 

{vt+vi) 



1.2 I-Q fOO 



0.1 -'0 J-00 



OQ 'J ~00 'J —00 



2 , _,^2'^ dp^{z,m) 
dm 



L' X 

•' if Lb'' SI 



f^'sS 'TisiOdCdzdv'^.dv'^Jv's^dv's 

»1.2 

--2+2/3, 



C+'^nsiOdC, 



(11) 



0.1 



where Ie appears only in the expression for m in 
dpL{z, m)/dm. 

As a check on the procedure, we eliminate z instead 
of m as an independent variable in Eq. (^). I.e., 
TO, Z, C, v;,, V'^a '"'sb' '"'si go to m,is, C, 
v'g^) with dz/dtE being the only term in the Jacobian 
determinant. Implicit differentiation of the square of 
Eq. (I) yields 



dz 

where 



^^(C-2z) _ tic f^dv^ ^ ^^dv, 

2z(C — z) mz{(, — z) \ '' dz ^ dz 



(12) 



dz 



dvi 
dz 



{-1 + Vm - WsJ/C- 



Some important assumptions are used in arriving at 
the above forms of these partial derivatives that will 
be pointed out below. From Eq. (H), we can write 
az^ -I- 6z -I- c = 0, with 

a = m + tl[{l-vsif +vl,]/C, 



b = -mC + 24K,(l 
c = tlC{vli+vl,), 



l) 



where Uv 



are transverse velocities rela- 



five to the observer. The solutions of the quadratic 
equation for z allow its expression in terms of the re- 
maining variables where it occurs in dz jdis and else- 
where in the integrand. The change of variables in 
the integrand of Eq. (^) is effected by replacing dz 
with [dz/ dtE)dtE- Integration over Ie fixes at the 
value selected by the 5 function, and the remaining 
six dimensional integral is 
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dz 
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0.1 



(13) 



With fixed tE, the solutions of the quadratic equation 
for z, the fractional distance of the lens to the center 
of the galaxy, for random choices of the independent 
variables within their physical ranges must be real 
with the added condition that z < (. Those sets of 
variables for which there are no real solutions for z 
are simply incompatible with the particular value of 
Ie- The two expressions for F, Eqs. ( |ll| ) and (p^, 
should yield identical time scale frequency distribu- 
tions. We shall demonstrate this below, but use the 
simpler form, Eq. (|l^) in subsequent calculations, 
wherein it is easy to introduce additional z and C de- 
pendence into the velocities. 

3. Models 

There is an arbitrary degree of complexity that 
one may incorporate into the galactic models that 
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specify the distributions of stars and their veloci- 
ties, and many of the parameters are not well con- 
strained by observations. The random velocities 
v'^^ f = {("^i ^1 ^7 ^) general with similar depen- 
dencies on parameters for v'^,^ ^ . The dependency of 
the velocity dispersion on stellar mass m is consistent 
with the increasing thickness of the spatial distribu- 
tion in the disk with decreasing mass {e.g. Mihalas 
and Binney, 1981). The best model of the galactic 
bar consistent with the COBE IR intensity distribu- 
tion (Zhao 1996) shows exphcitly the dependence of 
the dispersions on b and £. Part of the dependence 
of the velocities on z can be accounted for by sep- 
arating the galaxy into disk and bulge regions with 
different velocity dispersions in each region but with 
the dispersions being otherwise independent of z and 
C respectively for lenses and sources. This latter as- 
sumption was used explicitly in evaluating dvi^e/dz 
in Eq. (p^), where the z dependence of the velocities 
was ignored. We will relax this assumption in some 
of the examples below. 

The mean circular velocity of the gas clouds in 
the galactic disk appears to be a function of z or 
( (Clemens 1985), but the rather complex rotation 
curve is usually approximated as flat with all mean 
circular velocities equal to v^^r (Alvarez et al 1990). 
(Here we ignore any peculiar velocity of the LSR 
relative to the local mean circular velocity (Clemens 
1985).) The mean circular velocity of the stars is as- 
sumed to be the same as that of the gas clouds. A 
distribution of box-like orbits in the Zhao (1996) bar 
model matches the observed HI velocity distribution, 
and predicts the velocity dispersions that should be 
expected as a function of b and £ and presumably as 
a function of z or C inside the bulge. 

Microlensing should eventually prove a powerful 
constraint on the stellar mass function for the stel- 
lar spatial distribution that dominates the lens pop- 
ulation, which is apparently the central bar. This is 
especially so for the lower mass stars, since there is no 
other way of observing these stars. Generally, one ex- 
pects the mass function to depend on z (or C) in gen- 
eral, but that dependence is unknown. The frequency 
of event time scales is the microlensing constraint on 
the mass function, where the fraction of shorter time 
scales increases as the fraction of the mass density 
that is in the smallest stars or brown dwarfs is in- 
creased. However, increasing the velocity dispersions 
or otherwise increasing relative transverse velocities 
between lenses and sources and concentrating more 



of the lenses closer to the center of the galaxy both 
increase the fraction of short time scale events. So 
constraints on dispersions and spatial distributions of 
lenses and sources must be secure before definitive 
mass functions can be ascertained. 

Wc have adopted above the assumption of Kiraga 
and Paczyhski (1994) that the luminosity function of 
the sources is such that the fraction of stars with lu- 
minosities greater than some luminosity is pro- 
portional to . This appears to be a fairly good ap- 
proximation over the important range of luminosities 
in K band (Tiede, Frogel & Terndrup 1995), but less 
so in the visual band where the microlensing observa- 
tions are made. More complex luminosity functions 
may affect conclusions based on microlensing data 
but are not warranted by observational constraints 
at this time. Data provided in Tiede et al (1995) 
yield /3 « — 1 , which value is assumed in most of the 
examples herein. 

In principle, one could incorporate any of these 
complexities into the galactic model, but there is little 
motivation for doing very much of this until observa- 
tional constraints become more secure. Even without 
this security, one can determine the effect of various 
assumptions on the microlensing time scale rate and 
optical depth (Zhao and deZeeuw 1998) in anticipa- 
tion of future selection. For our purposes here, we 
shall adopt a galactic model for the integration of Eq. 
( prj ) that consists of a double exponential disk and a 
bar-shaped bulge. The disk model is that of Bahcall 
and Soneira, (1980). 



dp 
dm 



dm 



Pmo exp 
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300pc 



(14) 

where z' is the coordinate perpendicular to the plane 
of the galaxy, r is the radial coordinate in the plane 
of the galaxy, and is the scale length in the radial 
direction, which may be less than the 3.5 kpc used 
by Bahcall and Soneira (Zhao, Spergel & Rich 1995; 
Kent, Dame & Fasio 1995), and where /5„o is chosen 
such that /),„ = 0.05Mq/pc'^ at z' = and r = 8 kpc. 
Like Zhao et al (1996), we choose the triaxial bulge 
model of Zhao (1996) but keep the nucleus and do not 
terminate the bulge at 3.3 kpc. 
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where 
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(16) 



with g„ = 0.6, tr^ = 1.49 kpc, cr^ = 0.58 kpc and 
(Tj = 0.40 kpc. The coefficient p„ determines the mass 
of the bulge. The observations are consistent with 
the long axis of the bulge inclined about 13° to 20° 
relative to the line of sight to the galactic center with 
the near side of the bar lying in the first quadrant. 
There are large uncertainties in this inclination angle, 
however (Zhao 1998). 

In the applications, the origin of coordinates for 
both disk and bulge distributions will be transferred 
to the position of the observer with the following 
transformation. 



cos 9 — C cos b cos {I — 6) 
— sin 6* — C cos b sin {£ — 9 
C sin b 



(17) 



for a source in the bulge with xyz' being galacto- 
centric coordinates along the bulge principal axes and 
6 is the inclination of the bulge axis to the line be- 
tween the observer and the galactic center measured 
in the direction of increasing £. If a bulge star acts as 
a lens, ( is replaced by z in Eqs. ([l7|). For the disk 
distribution, the transformation is 



= 1 + z^cos^b- 2zcos6cos£ 



(18) 



Since b and £ are always small angles, cos b « cos £ fa 1 
is usually adequate approximation, and the exponent 
in Eq. jl^ ) can be replaced by (1 — s^Dg sin |6|/300pc) 

DgZ/Sa- 

In some of the examples below, we shall assume 
the circular velocity of the disk to be 210 km/sec 
independent of radial position r, with the random 
velocities superposed on this base. The bar can be 
assumed uniformly rotating at an angular velocity 
of 63.6km/sec/kpc, which is slightly higher than the 
60km/sec/kpc pattern speed chosen by Zhao (1996) 
but is arbitrarily chosen to match the disk rotation 
at 3.3 kpc. A bulge velocity dispersion will be su- 
perposed on this bar rotation. We shall also consider 
the case where the bar is not rotating, but neither 
of these assumptions is consistent with the results of 



Alvarez et al (1990) who find the flat rotation curve 
of the galaxy gas clouds (210 km/sec) to persist down 
to 2 kpc from the center with velocities exceeding 210 
km/sec at smaller galacto-centric distances. However, 
there are non circular streaming motions in the bar 
that could boost the apparent rotation inside 3.3 kpc 
(Zhao 1996). Finally, in most of our examples we shall 
consider the case where both circular velocities and 
velocity dispersions for lenses and sources are func- 
tions of where they are located in the galaxy. E.g., 
Circular velocity = 210 km/sec and velocity disper- 
sions are those appropriate to the disk for stars fur- 
ther than 3.3 kpc from the galactic center, but would 
be those appropriate to the bulge for smaller radial 
distances. 

Because there is no compelling argument specifying 
the mass function variation with galactic radius, we 
shall assume a single mass function applies through- 
out the galaxy. Most of the mass functions will be lim- 
ited to the range 0.08 < in < 2.0Mq. The upper limit 
is chosen because there are so few stars with larger 
mass, and the lower limit is chosen arbitrarily at the 
hydrogen burning limit. The consequences of using 
several different mass functions, sometimes extended 
into the brown dwarf region, will be demonstrated 
in examples that illustrate the sensitivity of the mi- 
crolensing technique to the mass function's eventual 
constraint. 

4. Monte Carlo integration 

For the choice of variables z, C, v^^^, Vli, Vsb, Vsi (Eq. 
(prj)) with no explicit z or C dependence for the dis- 
persions, the volume over which the integration vari- 
ables are distributed randomly in the Monte Carlo in- 
tegration is chosen to be (1.1)(1.2)(6'*) (JL^aLiasbCfsi 
so that the random velocities are truncated at ±3(Txb- 
For the examples where the dispersions and circular 
velocities are functions of z or C, all of the a's in the 
expression for the volume are replaced by the max- 
imum among them. Those random values of z > C 
are outside the volume over which the integration 
is taken. For the variables m, C, I'ib, t^if, Wst, '^sf (Eq. 
([!)), the volume is (1.92)(l.l)(64)a^,CT„CTs.CT«- We 
choose 0.1 < C ^ 1-2 since there are not likely to be 
sources closer than 0.8 kpc and faintness and blending 
will likely make sources beyond about 10 kpc unusable 
with the possible exception of the clump giants. The 
values of the six variables in the integrand are chosen 
randomly for each evaluation of the integrand, and 
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values of m in the first set and z in the second as 
functions of these variables and the fixed are de- 
termined by the solutions of Eq. (^). All points for 
which m is not within the specified range in the first 
form of F or where there are no real solutions for z 
(the particular choice of random variables not being 
consistent with the particular time scale) in the sec- 
ond or for z > C do not contribute to the integral. In 
the second form, the integral is evaluated separately 
for each of two viable values of z and the results added 
to get the total contribution to the integral for a par- 
ticular set of random variables. A sufficient number of 
points in the 6 dimensional volume are used in evalu- 
ating the Monte Carlo integral (typically 10^) to yield 
a relatively smooth frequency distribution as a func- 
tion of t^. 

5. Examples 

First we show that both means of calculating 
the time scale frequency distribution (Eqs. ( [Ill ) 
and (|l^)) yield the same result. For this purpose, 
we assume Vlo = v^sr = 210 km/sec, Vso = 0, 
(a^,,a„,crs.,c^«) = (16, 20, 124, 102) km/sec, {b,£) = 
(—4°, 4°), where the bulge velocity dispersions corre- 
spond to those given for these coordinates by Zhao 
(1996), bulge mass = 2.2 x 10^°Mq, local mass den- 
sity in the disk model = O.OSAfp/pc'^ with the Basu 
and Rana mass function (Eqs. (^^), and disk scale 
factors are 2.7 kpc and 300 pc. Fig. shows for val- 
ues of 6* = —20° and —65° that the two means of 
determining the distribution yield results that differ 
by amounts that arc at most comparable to the error 
estimates in the Monte Carlo integrals. This gives one 
some confidence in the correctness of the procedure. 
The first form of the distribution has a much sim- 
pler integrand and requires about a fourth as much 
computer time, so it will be used hereafter in all cal- 
culations. The 189 days and 12.6 million source stars 
used to scale the ordinate are the span of the MA- 
CHO (MAssive Compact Halo Objects) bulge obser- 
vations for the 1993 bulge season (Alcock et al 1996) 
and the number of sources they monitored. All of the 
time scale frequency distribution curves in the figures 
are such that /^^^^ F{tE)dtE is the number of events 
with Iei < Ie < tE2 for the 12.6 million source stars 
over the 189 day period. Alcock et al (1996) empiri- 
cally determine a sampling efficiency of detection as a 
function of event time scale, and they further reduce 
this efficiency by a factor of 0.75 for main sequence 
turnoff stars to account for the effects of blending. 
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Fig. 1. — Comparison of two Monte Carlo integra- 
tions with different variables eliminated from the in- 
tegrand. Zhao (1996) bar model with nucleus {M = 
2.2 X 10^°Mo). Double exponential disk model with 
scale lengths of 2.7 kpc and 300 pc and local mass 
density of O.O5M0/pc^. v^o = Wi,sr = 210 km/sec, 
vs„ = 0. (a^„(T„,cTs„as,) = (16,20,124,102) 
km/sec. The ordinate corresponds to 12.6 million 
source stars monitored over 189 days. 

There is no reduction of the sampling efficiency as- 
sumed for the clump giant sources. We use the recip- 
rocal of the overall efhciencies as the inferred number 
of events of a particular time scale that would have 
been detected had the overall efficiency of detection 
been 100%. E.g., about 7.41 4.4 day events would 
have been detected instead of one if 100% efhciency 
were possible. The MACHO detections as thereby 
enhanced are binned into 5 day intervals, and a his- 
togram of the distribution of the number of MACHO 
detections in the 189 day period of observation is plot- 
ted for comparison to the model distributions of time 
scale frequencies below. Only the 40 secure events for 
single stars are retained in the MACHO data set. 

For the remainder of our examples with the excep- 
tion of Fig. ^, we choose the velocities of source and 
lens to be determined by their position in the galaxy. 
For r > 3.3 kpc, we assume the stars have "disk" 
properties, with a flat rotation curve at 210 km/sec 
and velocity dispersions represented by a^t, = 16 
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km/sec and ct^^ = 20 km/sec in the directions of galac- 
tic latitude and longitude respectively. For r < 3.3 
kpc, the stars have "bulge" properties with circular 
velocities decreasing linearly from 210 km/sec at 3.3 
kpc to zero at the center of the galaxy, and with veloc- 
ity dispersions a^^ — <Ju — 110 km/sec. The dimen- 
sionless distances z and C are chosen randomly from 
their respective intervals and their values determine 
the velocity properties of the lens and source respec- 
tively. For example, both lens and source could lie 
in the bulge with high velocity dispersions for both. 
We shall also assume the Zhao (1996) bulge model 
described above {M — 2.2 x lO^^Mg) along with a 
double exponential disk model with scale lengths of 
2.7 kpc in the plane and 300 pc perpendicular to the 
plane for all cases below. We shall truncate the in- 
ner part of the otherwise singular nucleus in the Zhao 
model by requiring the stellar mass density to depend 
only on z' for values of r = + y'^ < 560 pc. This 
avoids the very large event rates and optical depths as 
lines of sight pass near the otherwise singular density 
at the galactic center. 

Our choice of the z and C dependence of the stellar 
velocities is somewhat arbitrary, although the rota- 
tion of the inner galaxy is close to the rigid body ro- 
tation of the bar pattern adopted by Zhao (1996), and 
the flat velocity curve of the disk is that determined 
by Clemens (1985) and Alvarez et al (1990). However, 
the latter two sets of authors find the flat rotation 
curve to apply for r > 2 kpc with gas cloud veloci- 
ties even increasing at smaller distances. These high 
central velocities may not be circular, however. The 
z and C dependence of the velocities used here is per- 
haps the simplest possibility that is semi-consistent 
with the observations, although more complicated 
(still unmotivated) dependences could be handled in 
principle. 

Six mass functions that will be considered are 
shown below. Proper overall coefficients are intro- 
duced in the integrals to yield the Bahcall and Soneira 
(1980) local density of O.O5M0/pc3 for the disk and 
a bulge mass of 2.2 x 10^°Mo. 

Basu & Rana 

(j){m) = 1.2588m"i 0.08 < m < 0.5327, 
= 0.3527m~3 0.5327 <m< 1.205, 



= 0.1401m- 



= 0.7294m" 
Modified Basu & Rana 



-6.83 



(m) = 0.1919m 



-2.5 



1.205 < TO < 2.0,(19) 



0.08 < TO < 0.5327, 



0.2861m 



-6.83 



Holtzman et 



(t){m) = 0.7795TO" 
= 0.5081TO" 

Modified Holtzman et al 

(j){m) = 0.2542to^2-^ 



0.5327 < TO < 1.205, 
1.205 < m < 2.0,(20) 



0.08 < TO < 0.7, 
0.7 < TO < 2.0, (21) 



0.08 < m < 2.0 



(22) 



Kroupa 



= 0.1038to-2-3^ 
= 0.6529to^-^ 
= 0.2674to-^-3^ 



0.075 < m < 0.35, 
0.35 < m < 0.6, 
0.6 < TO < 10.0, (23) 



Gould, Bahcall, & Flynn as extended by Mera et al 

0(to) = 4.55 X IO^^'to"'^ 0.08 < m < 0.11, 
= 0.823lTO""-^ 0.11 < m < 0.6, 
= 0.3371to-2-3^ 0.6 < to < 10.0. (24) 

The mass functions are such that m4>{'m)dm repre- 
sents the mass of stars/pc'^ in mass range dm about 
TO normalized such that the total mass density is 
IMq/pc^. The first mass function (Eq. [l9|) is that 



of Basu and Rana (1992) for local stars, and the sec- 
ond (Eq. (pO|)) is the same except the index on the 
low mass region is changed from -1 to -2.5. The third 
mass function Eq. ( pl| ) is determined by (Holtzman 
et al, 1998), and the fourth (Eq. ( p2[ ) is our modifica- 
tion, where we have changed the index from -1 to -2.2 
in the low mass region. The fifth and sixth mass func- 
tions (Eqs. (|2|) and (|4|) are those of Kroupa (1995) 
and Gould, Bahcall & Flynn (1997) as represented by 
Mera, et al (1998). The coefficients are those appro- 
priate to the lower mass cutoff of O.O8M0 except for 
the Kroupa mass function which is cut off at O.O75M0 
for later comparison with some results of Mera et al. 
We will extend the latter two mass functions into the 
brown dwarf region where the coefficients are altered 
to be consistent with the normalization of lA^o/pc'^. 
The modification of the Gould et al mass function by 
Mera et al consists of the top line in Eq. ( p^ ) for 
TO < 0.11. This modification is not supported by the 
Gould et al data, but we will keep it for the purposes 
of illustration. Also, the index of -0.6 in the second 



line of Eq (24) is reduced to about -1.0 if the unre- 
solved binaries are accounted for (A. Gould, private 
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communication, 1998). However, we shall maintain 
the -0.6 index again for illustration. The Kroupa and 
the Basu and Rana mass functions are based on mea- 
sured parallaxes of the nearby stars, the Gould et al 
mass function is derived from a luminosity function 
of more distant M stars (Mera et al 1998) and the 
Holtzman mass function is determined for bulge stars 
within Baade's window. The Basu and Rana mass 
function and more often its modified form will be used 
in most of the following examples 

Fig. II shows the effect of changing the luminosity 
index (3 from -0.5 to -2.0 superposed on the distri- 
bution obtained from the 1993 MACHO bulge data 
(Alcock et al 1996). (Although we shall see that 
the modeled distributions should be averaged over 
the MACHO field positions for comparison to the 
empirical distribution of time scale frequencies, the 
choice of a line of sight toward galactic coordinates 
(6,^) = (-4°,1°) (in Baade's window) yields distri- 
butions very close to this average shown in Fig. 11,) 



Fewer sources are visible as this index gets more neg- 
ative, and the expected decrease in the number of 
events observed as the index varies from -0.5 to -2.0 is 
displayed. The modified Basu and Rana mass func- 
tion is used, and the curve for (3 — —\ is the same 
as that shown as the "complete model" in Fig. ^. 
The other model parameters are specified in the fig- 
ure caption. Fig. ^ shows the distribution of time 
scale frequencies for several bar orientations using the 
Basu and Rana mass function. The number of events 
decreases as the magnitude of the inclination of the 
bulge increases as expected, since fewer stars are lo- 
cated along the line of sight as the long axis of the 
bulge is pointed further away from the observer. Since 
the bulge distribution is so dominant over that of the 
disk, the time scale frequency distribution is almost 
proportional to the mass of the bulge. Hence, the 
curve for 9 = —13° is in some ways a closer match 
to the MACHO data if the bulge mass is increased 
by a factor of 1.5. Zhao et al (1996) also pointed out 
the sensitivity of the event rate to the bulge mass, 
but notice here that simply increasing the bulge mass 
to match the peak in the distribution produces too 
many events in the 8 to 40 day range. None of the 
model distributions seem to produce as many long 
time scale events as the MACHO group observed. To 
show how a different velocity model can affect the 
model time scale frequencies, we repeat the exercise 
of Fig. |[ where all parameters are identical to those 
used to construct that figure except the velocities. In- 
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Fig. 2. — Effect of changing the source cumulative 
luminosity function index (3. The decrease in the 
number of events as the magnitude of the index in- 
creases is as expected since fewer sources are visible. 
Zhao (1996) bar and nucleus model with bulge mass 
~ 2.2 X IQ^^Mq. Double exponential disk model 
with scale lengths of 2.7 and 0.3 kpc in the disk 
and perpendicular to the disk with local value of 
Pm = O.O5M0/pc^. Wio, Wso = -210 km/sec for z, C < 
0.5875, = 2.424(2, C - 1) for z, C > 0.5875. Modi- 
fied Basu and Rana mass function (0.08 < m < 2.0). 
{(Jab,aM,(Jii,aw) = (16, 20,110, 110 km/sec). The or- 
dinate corresponds to 12.6 million source stars moni- 
tored for 189 days. 



stead of assigning circular velocities and velocity dis- 
persions as a function of position in the galaxy, we 
here assume Wi,o = 210 km/sec for all z, Vga — and 
(Tj^b = a^^f = Usb = o'st = 110 km/sec for all z and 
Although this dispersion is not realistic, it might 
be a reasonable approximation if most of the lenses 
were in fact in the bulge. On the other hand, keep- 
ing the source circular velocities zero while assuming 
that all of the lenses maintain the flat velocity curve 
of the disk already introduces a significant relative 
transverse velocity between source and lens before 
the dispersions are applied. Then the large disper- 
sions superposed onto this velocity difference result 
in a large increase in the relative number of short 
time scale events and in the total number of events. 
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Fig. 3. — Effect of changing the inclination of the bar 
axis to the Hne of sight. The model is the same as 
that in Fig ^ except the unmodified Basu and Rana 
mass function is used. 

The results shown in Fig. ^ show this large increase 
in the time scale frequency distribution for all of the 
bulge orientations. Even with the use of the Basu 
and Rana mass function, the MACHO data is now 
well matched by the distribution for 9 = —13° with- 
out increasing the bulge mass, but this match is an 
artifact of an unrealistic velocity model. We return 
now to the velocity model of Fig. ^ where both circu- 
lar velocities and the dispersions depend on z and Q 
along with the other parameters assumed in the con- 
struction of that figure to show explicitly the depen- 
dence of the model time scale frequency distribution 
on velocity dispersion. Fig. ^ displays the distribu- 
tions for several choices of velocity dispersions, where 
more long time scale events are produced but with 
fewer total events, and the peak in the distribution 
moves to longer time scales as the velocity dispersion 
is decreased. Notice, however, that increasing the ve- 
locity dispersion of the disk part of the galaxy to the 
values appropriate to the bulge results in only a mod- 
est increase in the distribution over that where only 
the bulge has the high dispersion. This is because 
in the latter case most of the lenses already have a 
large dispersion, since relatively few lenses are in the 
disk distribution. The sensitivity of the microlensing 
time scale frequency distribution to the stellar mass 



Fig. 4. — Effect of changing the velocity model. 
Galactic model is the same as Fig|| except v^q — 210, 
Vsa — km/sec, = = <Jsb = cfsi = 110 km/sec. 
Short time scale event rate and total number of events 
is increased significantly. 

function of the lenses is demonstrated in Figs. ^ and 
0, where consequences of assuming the several mass 
functions given in Eqs. ( ]T9| ) to ( p^ ) for various ex- 
tensions into the brown dwarf region are compared 
with the MACHO data. In Fig. | we show the time 
scale frequency distributions with 9 = —13° for three 
Kroupa mass functions with successively larger excur- 
sions into the brown dwarf region as indicated by the 
minimum mass. The indices shown in the figure ap- 
ply only for the lowest mass range in the definitions 
of the mass functions. Except for the mass functions, 
all model parameters are those of Fig. ^ with (3 = —1. 
Fig. shows the distributions of the Gould et al 
mass functions, where the same model assumptions 
apply as in Fig. H. First notice the degree of sensi- 
tivity of the model time scale frequency distributions 
to the mass functions, the only variation in the mod- 
els. Herein lies the eventual power of the microlensing 
technique to constrain the galactic mass function far 
from the Sun. The Kroupa mass function extended 
down to m = 0.05 yields a distribution that is a fair 
match to the MACHO data, although it is somewhat 
"thin" in the region around the peak. Improving the 
fit on the longer time scale side of the peak by in- 
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Fig. 5. — Effect of changing the velocity dispersions. 
The model is the same as that in Fig. || except for 
the changed velocity dispersions, the fixed value of /3 
and the use of the unmodified Basu and Rana mass 
function. 

creasing the galactic bulge mass makes the fit worse 
for near 33 days. A similar problem is faced by the 
Basu and Rana mass function for the bulge mass in- 
creased by a factor of 1.5 in Fig. ^, where there are too 
many events in the 8 to 40 day range of tj^. Perhaps 
the best fit to the MACHO data is for the modified 
Basu and Rana mass function in Fig. || for /? = — 1 
and = —13° even though the peak is somewhat too 
high. This fit could even be improved somewhat by 
decreasing the mass of the bulge slightly (this would 
lower the peak), while making the -2.5 index in the M 
star region slightly less negative (this would "fatten" 
the curve to bring it back into coincidence with the 
Macho data. It is noteworthy that the best fit to the 
MACHO data is obtained with only hydrogen burn- 
ing stars, albeit with the M star region enhanced with 
an index of -2.5. 

Mera et al (1998), with a different galactic model, 
find a time scale frequency distribution that is consis- 
tent with the MACHO data using the Kroupa mass 
function with lower bound O.OTSM© — slightly into the 
brown dwarf region. This same mass function yields 
a poor match to the MACHO data with the model 
used here. This illustrates the model and perhaps 



Fig. 6. — Comparison of the time scale frequency dis- 
tributions obtained from the Kroupa mass functions 
extended into the Brown dwarf region with that from 
the MACHO 1993 bulge data. All model parameters 
except the mass function are those of Fig. ||. 

procedure dependence of conclusions about the mass 
functions when so may crucial galactic parameters are 
so poorly constrained by observations. 

Extending the Kroupa mass functions further into 
the brown dwarf region to enhance the time scale fre- 
quency distribution is only partially successful. If one 
believes the statistics of small numbers characterizing 
the MACHO data, increasing the fraction of the mass 
distribution in the brown dwarfs makes the distribu- 
tion too narrow and shifts the peak too far toward 
shorter time scale events. The distributions so de- 
rived have (slightly) too few events of moderate time 
scales. The Gould et al mass functions, (Fig. ^ are 
even worse matches to the MACHO data as also found 
by Mera et al for their model. This mass function has 
many fewer stars between 0.11 and O.35M0, and for 
m > 0.08 there is only a short range where the mass 
function appears to turn up steeply in the Mera et al 
extension. The severe lack of small mass stars leads to 
the shift in the peak in the distribution to longer time 
scales, more moderately long period events, and a 
rather modest total number of events. Extending the 
very steep mass function into the brown dwarf region 
(rn > 0.04) fails worse than the same exercise with 
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the Kroupa mass function in matching the MACHO 
data since the peak of the distribution is moved even 
further toward very short time scales (< 4.5 days), ex- 
ceeding the MACHO numbers there while failing by 
a factor of at least two to produce a sufficient number 
of events of time scales between 5 and 20 days. Since 
the modified Basu and Rana mass function gives the 
best fit to the MACHO data, we shall use this mass 
function in subsequent examples, except in Fig. |ll| 
where the modified Holtzman et al mass function is 
shown to also give a reasonably good fit to the MA- 
CHO data. Fig. ^ shows the relative importance of 
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Fig. 7. — Same as Fig. ^ except Gould et al mass 
functions are used instead of the Kroupa mass func- 
tions. Note, however, that the extension of the Gould 
et al mass function into the brown dwarf region by 
Mera et al (1998) is completely speculative. 

various contributions to the model mass distribution 
in determining a time scale frequency distribution for 
9 = —13°. Removing the disk stars has relatively 
little effect on the distribution confirming the suspi- 
cion that most of the lenses are in fact located in the 
bulge. Consistently, the number of events is drasti- 
cally reduced when only the disk is present, although 
the peak in the distribution is not changed much be- 
cause the velocities still have the z and C dependence 
of Figs. ^ and ^. Removing the nucleus decreases the 
number of events by a factor of about 2.5. The nu- 
cleus contributes about 21% of the total bulge mass. 



but has such a large effect because the line of sight 
to {b,£) = (-4°, 1°) for e = -13° encounters a com- 
parable number of nucleus and bar stars. We should 
point out here that Zhao et al (1996) are able to get 
a good match to the MACHO data with a compara- 
ble mass function and velocity dispersion as that used 
here without the nucleus, whereas we require the nu- 
cleus be kept for a good match. Like the Mera et al 
(1998) match to the MACHO data with the Kroupa 
mass function which failed for our model, this again 
illustrates the model and procedure dependence of 
estimated time scale frequencies that is independent 
of and thereby weakens the constraints on the mass 
function. The Zhao model is also very sensitive to the 
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Fig. 8. — Relative contributions of several model con- 
stituents to the model time scale frequency distribu- 
tions, where the model is that of Fig. |^. 

galactic coordinates of the line of sight. This is shown 
in Figs. PI and [l0|, which show the time scale frequency 
distributions along two perpendicular traces through 
Baade's window whose galactic latitude (—4°) bisects 
the more extended region where the MACHO dat was 
obtained. The smaller scale lengths assumed in the z 
direction for both the bar and nucleus results in the 
more rapid fall off in events at higher galactic lati- 
tudes than for higher galactic longitudes. Recall that 
we have eliminated the singular nature of the nucleus 
by truncating the nucleus density at the galactic cen- 
ter by assuming that density uniform for constant z' 
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inside 560 pc. The nucleus thereby loses about 36% 
of its mass but the whole bulge loses only 7%. If 
the nucleus is not truncated, both the rate of events 
and the optical depth become very large as the line of 
sight approaches the center of the galaxy. Zhao et al 
(1996) avoid this problem by eliminating the nucleus 
altogether, but we require the nucleus stars to match 
the MACHO data as discussed above. 
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Fig. 9. — Dependence of the time scale frequency dis- 
tribution on the galactic latitude of the line of sight 
with the longitude fixed at 1°. The singularity in the 
nucleus density distribution is removed as discussed 
in the text. Model is that of Fig. ||. 

The large dependence of the event rate on coor- 
dinates of the line of sight means that microlensing 
will also eventually prove a powerful constraint on 
bar model star distributions as well as mass functions 
if other parameters such as the velocity dispersions 
can be adequately defined. As the MACHO fields 
are spread over several degrees near Baade's window, 
(Fig. 1 of Alcock et al (1996)), it is appropriate that 
the model event frequency distributions be averaged 
over the MACHO fields for comparison with their re- 
sults. This average is shown in Fig. |ll| for Basu and 
Rana and Holtzman et al mass functions and their 
modifications by the more negative indices for low 
mass stars. Both modified mass functions yield good 
agreement with the data with the best match prob- 
ably obtained with an index between -2.2 and -2.5. 
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Fig. 10. — Same as Fig. || but longitude is varied at 
constant —4° latitude 

However, given all of the uncertainties in the model 
parameters discussed above and the simplicity of the 
model itself, this agreement must be regarded as for- 
tuitous and should not be regarded as a definitive 
constraint on the mass function. Note, however, that 
there is as yet no need to assume that there is any sig- 
nificant fraction of the bulge mass in brown dwarfs — a 
conclusion also reached by Zhao et al (1996). On the 
other hand, we do need more M stars in the mass 
functions than obtained from the star counts in or- 
der to match the MACHO data. There may be a 
not unreasonable galaxy model that will yield a suf- 
ficient number of short time scale events with a m^^ 
variation of the mass function in the M star region, 
but the necessary changes in the parameters of the 
model adopted here seem too extreme. But the un- 
certainties in the most relevant mass function of the 
galactic bulge stars obtained by Holtzman et al, the 
paucity of the microlensing data and modest observa- 
tional constraints on the galactic model itself do not 
yet warrant a search for a drastically different galactic 
model. We now turn to a determination of the opti- 
cal depth as a function of the line of sight. We shall 
see that a similar average over the MACHO lines of 
sights for our best matching model also gives reason- 
able agreement with the MACHO estimates. We shall 
eventually conclude that the Zhao model distribution 
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Fig. 11. — Model time scale frequency distributions 
averaged over the lines of sight to the 24 MACHO 
fields for modified and unmodified Basu and Rana 
mass functions and modified and unmodified Holtz- 
man et al mass functions. The variation of the mass 
functions in the M star region is indicated in paren- 
theses. The model is that of Fig. || except for the 
coordinates of the lines of sight and the explicit mass 
functions. 



of stars with our assumptions about the mass function 
and the velocity distribution yields a good match to 
both the estimated optical depth and the time scale 
frequency distribution of the 1993 MACHO data, but 
that the meagerness of this data and the uncertainties 
within the galactic models means that the microlcns- 
ing constraints on both the mass function and the 
distributions of lens and source stars are very weak 
indeed. 

6. Optical depth 

If we use the definition of a microlensing event as 
a source at distance Dqs passing within the Einstein 
ring radius of a lens at distance I?oli the cross 
section of the lens for such an event is just 7ri?|, 
where RE{m, Dol, Dos) is given by Eq. (|). For 
such lenses within a slab of thickness dDoL at Dql, 
the probability dP that a ray from a source will pass 
through an Einstein ring on its way to the observer 
is thus TrR^{dnL{m, DoL)/d'm)dmdDoL, where Ul is 



the number density of lenses and where we have se- 
lected only the lenses with masses within dm of m. 
If / is the intensity of rays from the source that have 
not passed through an Einstein ring in traveling dis- 
tance Dqs — Dq^ on their way to the observer, then 

oi from lenses 
dP or 



the fractional change in traversing dD 
with masses within dm oi m is dl / 1 



dl 

T 



AirGm Dol{Dos - Dol) 1 dpL{m,DoL) 
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(25) 

where is the stellar mass density and Eq. has 
been used. Integration of Eq. ( p5| ) leads to the optical 
depth for microlensing to distance I?os of 
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where cancellation of the masses in Eq. ( |25| ) allows 
immediate integration over the lens mass function. 
Finally we introduce the dimensionless variables used 
earlier and average the optical depth over the dis- 
tribution of visible sources along the line of sight to 
obtain the result of Kiraga and Paczynski (1994) 
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where ^^(C) is the source number density and where 
the coefficient is 3.86 x 10~^ if p^ is expressed in 
Mq/pc^. We shall omit the brackets around r here- 
inafter, but always understand that it is averaged over 
the source distribution. 

The integral in the numerator of Eq. ( p7[ ) is eval- 
uated with a Monte Carlo technique like that used 
for the time scale frequency distribution, whereas the 
denominator is evaluated with a standard numerical 
technique. Our examples of optical depth will be 
mostly limited to the model of Fig. ^, as that model 
which best matches the time scale frequency distribu- 
tion averaged over the MACHO lines of sight. This 
model is the Zhao bar and truncated nucleus with to- 
tal mass of 2.2 x and the double exponential 
disk with scale lengths of 0.3 and 2.7 kpc, but the op- 
tical depth is independent of the mass function. The 
optical depth is quite sensitive to the coordinates of 
the line of sight for this model as indicated in Fig. 
|l2|, where we show r as a function of latitude and 
longitude around Baade's window. The dots in Fig. 
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mark the coordinates where we have evaluated the 
time scale frequency distribution in Figs. |^ and |l^, 
and they display the expected correlation between the 
optical depth and number of events. If we reduce the 
radius inside of which the mass density of the nu- 
cleus is constant, the optical depth, like the time scale 
frequency distribution, is increased for latitudes and 
longitudes less than about 4°, but is little affected 



for larger coordinates. In Fig. 13 optical depths for 
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Fig. 12. — Dependence of optical depth on direction 
of the line of sight for the Zhao bar and truncated nu- 
cleus (M = 2.2 X IO^'^Mq) with double exponential 
disk with scale lengths of 0.3 and 2.7 kpc and local 
density of O.O5M0/pc"^. The dots correspond to the 
directions for which the time scale frequency distri- 
bution was determined in Figs. and HG. 



with the Basu and Rana mass function by a factor of 
1.5 to match the peak in the empirical distribution. 
Although this leads to too many mid duration events, 
it could be regarded as a better fit. This multiplica- 
tion simulates a similar increase in the bulge mass. In 
Fig. |l^, we see that the optical depths for the MA- 
CHO field coordinates are more compatible with the 
empirical estimate for the higher mass bulge. How- 
ever, in addition to the poorer match of the Basu 
and Rana mass function to the moderate duration 
events obtained by MACHO in Fig. ||, the higher 
mass bulge may be less compatible with the stellar 
kinematics as determined by Zhao (1996). The model 
optical depth averaged over the MACHO field direc- 
tions is 1.54 X 10"*^ for a bulge mass of 2.2 x W°Mq 
and 2.14 x 10"'^ for 3.3 x 101°Mq compared with 
(2.4 ± 0^ X 10-s for the MACHO empirical value. 
In Fig. |j we show the effect of different positions of 



o 

X 3 



o 



- 1 ' ' ' 


1 1 1 1 1 1 

e=-i3" 






/? = -! 




2a 












- 1(7 













MACHO T 




















* t * ♦ 










- . M = 2,2xl0 






- . M = 3.3xl0 




°M0 





2 4 6 

Galactic latitude b| (degrees) 



the MACHO line of sight directions are distributed 
across horizontal lines designating the estimated MA- 
CHO empirical optical depth averaged over all of the 
24 fields and 1 and 2ct error estimates. The points for 
a bulge mass of 2.2 x 10^°Mq, that which gave such 
a good match for the time scale frequency distribu- 
tion, fall somewhat low with about half falling within 
the 2(7 lines. In some sense, this can be regarded as 
not unreasonable agreement given the uncertainties 
in the galactic model and the sparseness of the data. 
Recall in Fig. || that we multiplied the time scale 
frequency distribution obtained for the Zhao model 
bulge together with our standard velocity model but 



Fig. 13. — Comparison of model optical depths for 
the MACHO fields with the empirically determined 
average value. The model is the same as that in Fig. 
[l2| except for the specific bulge masses indicated for 
the two sets of points. 

the bar axis along with luminosity function index (3 
for the direction (6,^) = (-4°, 1°). The optical depth 
shows the expected decrease as the magnitude of the 
index increases — location of fewer visible sources at 
the larger distances pushes the mean position of the 
sources closer to the observer with less optical depth 
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resulting. The optical depth also shows the expected 
decrease as the bar is rotated further from the line of 
sight. We show the optical depth for the near side of 
the bar in the fourth quadrant {9 — +13°) to demon- 
strate the expected decrease when the near side of 
the bar is not in the same quadrant as the line of 
sight. The dots indicate values of /3 and 9 for which 
we determined the time scale frequency distribution 
in Fig. 1^ where the optical depths and the number of 
events are again appropriately correlated. As a com- 
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Fig. 14. — Variation of optical depth with inclina- 
tion of the bar axis and with the source luminosity 
function index. The model is the same as in Fig. |l^ 
except for the indicated variations in 9 and /3. The 
dots correspond to the values of 9 and /? for which 
time scale frequencies were determined in Fig. ^. 

panion to Fig. || showing the relative contributions of 
various parts of the galactic model to the time scale 
frequency distribution, we show in Fig. |l^ the effect 
of the same dissection of the model on optical depth. 
Again, the disk contributes relatively little to the op- 
tical depth — like its meager contribution to the time 
scale frequency distributions. The nucleus becomes 
more important when the line of sight passes very 
close to the galactic center. The optical depth for the 
complete model would turn up sharply at small galac- 
tic latitudes in Fig. |l^ without the truncation of the 
nucleus. 
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Fig. 15. — Contributions of various parts of the model 
of Fig. ^to the optical depth as a function of galactic 
latitude. 

7. Discussion 

We have developed a procedure for determining 
model time scale frequency distributions for microlens- 
ing events toward the center of the galaxy that can 
be easily modified for improvements in galaxy models 
as they become available. Two methods of evaluat- 
ing the integral used to accumulate all of the com- 
binations of variables leading to the same time scale 
give the same distribution within the errors of the 
Monte Carlo integrations (Fig [|) . This robustness of 
the model distribution gives one some confidence that 
the computational procedure is sound. That changes 
in the distributions as various parameters are varied 
in a series of examples were as expected increases this 
confidence. 

The complete variable dependence of the event 
time scale is maintained within the integral used to 
accumulate all the combinations of variables leading 
to the same time scale. This contrasts to some work 
in the literature where event rates are determined for 
all of the lenses being IMq with a later convolution 
over the mass function (e.g. Kiraga and Paczynski 
1994; Zhao et al 1996). However, it seems implicitly 
assumed in this latter procedure that galactic mod- 
els are secure enough that the mass function can be 
separated out as the major unknown and mcaning- 
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fully constrained by substituting various forms into 
the convolution. There would also appear to be some 
danger of introducing some artifacts into the compu- 
tation by such piecewise averaging. Within the Zhao 
(1996) bulge model (bar plus nucleus) and the dou- 
ble exponential Bahcall and Soneira disk model that 
we have adopted, the time scale frequency distribu- 
tions are very sensitive to the total bulge mass, the 
circular velocity model, velocity dispersions and the 
direction of the line of sight as well as the mass func- 
tion. Other bulge and disk distributions of stars will 
introduce still other variations. So until the prop- 
erties of the galaxy are better constrained, it seems 
premature to assign the highest priority to the mass 
function in matching the model time scale frequency 
distributions to the observations. 

With our nominal adornments to the Zhao bulge 
model, we obtain a good match to the MACHO time 
scale frequency distribution with an optical depth 
that is only slightly too low (Figs. |ll| and |l^ ). But 
we do this only by changing the index of either the 
Basu and Rana or Holtzman mass function from -1 to 
-2.5 and -2.2 respectively in the M star region (Fig. 
pi] ). The Kroupa mass function has a slope of -2.35 for 
the latest M stars, but the Gould et al mass function 
is even flatter than that of Basu and Rana (-0.6 vs 
-1.0), with the steep extension into the brown dwarf 
region by Mera et al (1998) being completely specu- 
lative. From Figure |ll| one is led to an index some- 
where between -2.2 and -2.5 for an even better fit to 
the MACHO data within the overall model. But the 
Holtzman et al mass function for bulge stars in the M 
star region (~ m~^) where most of the lenses reside is 
like that of Basu and Rana and also Gould et al if un- 
resolved binaries are accounted for in the latter. So 
the only observational support for the steeper mass 
function in the M star region is consistency with a 
meager MACHO data set. Still, microlensing may be 
detecting M stars that are not accounted for in the 
star counts because of difficulties in estimating the 
completeness of the optical surveys (Holtzman, pri- 
vate communication, 1998), or the meagerness of the 
microlensing data set may make the empirical event 
rates not entirely representative. 

With a different velocity model, or with other than 
a Zhao bulge distribution of stars, our match may be- 
come less striking. Still, it is significant that we found 
a model consistent with the microlensing data with- 
out resorting to a large population of brown dwarfs 
as advocated by several authors {e.g. Han and Gould 



(1996); Han 1998). In fact, we find that extending the 
Kroupa and Gould et al mass functions into the brown 
dwarf region produced too many very short time scale 
events and an insufhcient number of moderate time 
scale events (Figs. ^ and although at least the 
Kroupa mass function including brown dwarfs might 
have been a better fit if we had just decreased the 
velocity dispersions or changed the bulge mass. 

Our use of the Kiraga and Paczyhski luminosity 
function index /3 as a measure of source visibility is 
useful for verifying consistent behavior of the distribu- 
tions derived by this procedure (Fig. ^), but it is actu- 
ally not a good representation. First, the luminosity 
function is not so simply behaved in the visible re- 
gion, and second, blending is a much more important 
limitation on source visibility than apparent bright- 
ness. The representation might be useful in account- 
ing for source visibility in highly obscured regions, but 
there are likely to be too few visible sources in these 
circumstances to motivate an extensive microlensing 
search. Penetration of these regions by using infrared 
wavelengths would return us to blending as the most 
important limitation. Although there is a need for a 
more realistic model of source visibility, it is probably 
less important than other uncertainties in the galactic 
model. 

The change in the distribution function as we al- 
tered the inclination 6 of the bar axis to the observer- 
galactic center line also showed consistent behavior 
for the model (Fig ^) as did the optical depth (Fig. 
[T^ ). Although this inclination is uncertain to some 
extent, which fact motivated our investigation of its 
effect, changing the inclination also requires chang- 
ing the shape of the bar to remain consistent with 
the COBE observations (Zhao 1998). We left the bar 
shape unaltered during this exercise, but other uncer- 
tainties preclude a more complicated model. 

It not clear just how well our velocity model rep- 
resents the real world, but a more complicated one is 
not warranted now. Certainly, the velocity dispersion 
must have a more complicated z and C dependence 
than the one that we used mostly {e.g. Figs. || and ||), 
and the flat circular velocity of the disk and ramped 
circular velocity of the bulge can only be rather crude 
approximations. Still, it is useful to realize the ex- 
treme variation in the calculated time scale frequency 
distributions simply by changing the velocity disper- 
sions within the model (Fig. ||) or by a change in the 
model itself (Fig. A velocity model, both circular 
velocities and dispersions or even more complicated 
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models accounting for the box like orbits in the Zhao 
(1996) bulge model, that does not represent the real 
galaxy very well may lead to, say, the adoption of a 
mass function in matching microlensing data that is 
far from reality. 

We have so far in this discussion concentrated on 
the uncertainties in deducing galactic properties from 
microlensing time scale frequency distributions and 
optical depths toward the galactic center. However, 
the sensitivity of such microlensing data to changes in 
the model illustrated by our series of examples shows 
just how powerful the technique will eventually prove 
to be in the study of galactic structure. Figs. ^, 
p^ , and |l3| show the very rapid variation in both 
optical depth and time scale frequency distributions 
with relatively small changes in the line of sight di- 
rection. If such variations do not emerge with the 
collection and analysis of many more events to ob- 
tain the necessary spatial resolution, one would have 
to drastically revise current ideas about the nature 
of the bar. Much better spatial resolution may al- 
ready be possible from the more than 200 additional 
events that have been monitored since 1993, but that 
remain unanalyzed and unpublished. Recall also that 
all ways to improve the fit to the data are not equiva- 
lent. In Fig. H, the Basu and Rana mass function led 
to a time scale frequency distribution that produced 
far too few short time scale events for 6 = —13°, one 
of the preferred inclination angles for the bar axis. 
We increased the number of these events for this in- 
clination by increasing the galactic mass by a factor 
of 1.5, an increase that might be tolerated by other 
observational constraints. If the MACHO distribu- 
tion is close to the correct one, the number of events 
predicted between about 15 and 40 days is now too 
large. A better way to match the data was to increase 
the magnitude of the index in the M star range of the 
Basu and Rana or Holtzman et al mass functions as 
shown in Fig. |ll|. 

So microlensing will place constraints on the galac- 
tic distribution of stars, their distributions of veloc- 
ities and the mass function, but other more conven- 
tional observations must also be accumulated to in- 
dependently constrain the parameters for a gradual 
convergence to the true properties of the galaxy. Cer- 
tainly with only 40 data points, very uncertain veloci- 
ties and still uncertain star distribution in spite of the 
success of the Zhao (1996) model (We have not even 
discussed deviations from axial, radial or bar sym- 
metry.), it is premature to assert any constraints on 



galactic properties in general from the microlensing 
data (except for the existence of the bar) and in par- 
ticular on the mass function. We found a mass func- 
tion that works along with our other assumptions, but 
we are unable to assert that it is the right one. 

8. Conclusions 

We end by enumerating our conclusions. 

1. The procedure developed here appears to be com- 
putationally sound (Fig. ^ and later consistencies) 
and provides a convenient framework for rapidly de- 
termining the predictions for microlensing surveys of 
the galaxy. Position dependence of various galactic 
properties such as stellar velocities is relatively easily 
incorporated when observations justify such detail. 

2. The sensitivity of the predictions of our overall 
model to changes in various galactic parameters used 
in the model, means that microlensing will be a major 
player in the set of many types of observations that 
will ultimately constrain the structure, content and 
dynamics of the galaxy. 

3. There is no motivation for invoking a large pop- 
ulation of brown dwarfs to account for the number 
of short time scale events observed by the MACHO 
group. Brown dwarfs might be accommodated by re- 
ducing the velocity dispersion, for example, but there 
appears to be a sharp fall off in the number of brown 
dwarfs as companions to local solar type stars (Mazeh, 
Goldberg & Latham 1998), and there is no appar- 
ent theoretical or observational reason for expecting 
many brown dwarfs in the bulge. 

4. The sensitivity of the microlensing event rates and 
optical depth to the direction of the line of sight for 
bar models means it will be a sensitive probe of the 
star distribution in such models. 

5. The disk stars contribute to a minor fraction of 
the microlensing events (Figs. || and p^ ). 

6. We obtain a better fit to the MACHO data by in- 
cluding the nucleus in the Zhao (1996) bulge model, 
albeit we truncated the singular peak in the density to 
avoid excessively large optical depths and event rates 
for lines of sight near the galactic center. This trunca- 
tion is probably academic in any case as obscuration 
increases drastically for lines of sight very close to the 
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galactic center and no data have been taken there. 

7. Our model time scale frequency distributions dance 
around the empirical distribution from the MACHO 
observations as various parameters are varied, but we 
are able to get a very good match to this data from the 
distributions from our nominal model averaged over 
the MACHO field directions (Fig. About half of 
the optical depths for these directions fall within the 
2ct error bars of the empirical average optical depth 
(Fig. |l^), and the mean optical depth is about 2a 
away from the MACHO mean. It is important to re- 
member that although no brown dwarfs are necessary 
to match the MACHO data, the number of M stars 
had to be increased rather drastically above those ob- 
tained from recent star counts both locally as in the 
Basu and Rana or Gould et al mass functions and in 
the galactic bulge itself from the Holtzman et al mass 
function. 

8. For our particular model, the addition of brown 
dwarfs to the mass function seems to be a detriment 
rather than an asset in matching the MACHO data 
(Figs. |and|7l). 

9. In spite of the success of our nominal model in 
matching the MACHO data, uncertainties in all of 
the galactic parameters and the meagerness of the 
published data preclude definitive microlensing con- 
straints on any of them — including the mass func- 
tion. This frustration should be diminished as the 
microlensing data set grows. In fact, additional exist- 
ing data may already yield sufficient spatial resolution 
to constrain the bar shape and orientation — if the set 
were analyzed and published. 
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